What can we learn from the Interacting Boson 
Model in the limit of large boson numbers? 

Dennis Bonatsos'^, E. A. McCutchan^, R. F. Casten*^ 

qq", ^Institute of Nuclear Physics, N.C.S.R. "Demokritos", GR-15310 Aghia 

^D ' Paraskevi, Attiki, Greece 

O ! u 

p^ , Physics Division, Argonne National Laboratory, Argonne, Illinois 60439, USA 



00 
(N 



^ 

a 



oo 

o 



^ Wright Nuclear Structure Laboratory, Yale University, New Haven, Connecticut 

06520-8124, USA 



Abstract 



O ' Over the years, studies of collective properties of medium and heavy mass nuclei 



in the framework of the Interacting Boson Approximation (IBA) model have focused 
on finite boson numbers, corresponding to valence nucleon pairs in specific nuclei. 
Attention to large boson numbers has been motivated by the study of shape/phase 
^ ' transitions from one limiting symmetry of IBA to another, which become sharper in 

^«0 , the large boson number limit, revealing in parallel regularities previously unnoticed, 

although they survive to a large extent for finite boson numbers as well. Several 
Tlj- ■ of these regularities will be discussed. It will be shown that in all of the three 

[->^ . limiting symmetries of the IBA [U(5), SU(3), and 0(6)], energies of 0"'" states grow 

^D \ linearly with their ordinal number. Furthermore, it will be proved that the narrow 

transition region separating the symmetry triangle of the IBA into a spherical and 
a deformed region is described quite well by the degeneracies £'(0^) = E{Q'\^), 
£(0+) = ^(10|), ^(0|) = -E(14+), the energy ratio £;(6^)/£;(0^) turning out to 
^ ' be a simple, empirical, easy-to-measure effective order parameter, distinguishing 

j^ ■ between first- and second-order transitions. The energies of O"*" states near the point 

of the first order shape/phase transition between U(5) and SU(3) will be shown 
to grow as n(n-|-3), where n is their ordinal number, in agreement with the rule 
dictated by the relevant critical point symmetries studied in the framework of special 
solutions of the Bohr Hamiltonian. The underlying dynamical and quasi-dynamical 
symmetries are also discussed. 



Collective phenomena in atomic nuclei are described in terms of two comple- 
mentary models, the algebraic Interacting Boson Approximation (IBA) model 
[1], and the geometrical collective model [2|3] . In the former, s and d bosons 
(bosons of angular momentum and 2 respectively) are used, while in the 
latter the collective variables (3 (the ellipsoidal deformation) and 7 (a mea- 
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Fig. 1. IB A symmetry triangle with the three dynamical symmetries. The critical 
point models E(5) and X(5) are placed close to the phase transition region (slanted 
lines). The solid curve indicates the Alhassid-Whelan arc of regularity. Adopted 
from Ref. [5j. 

sure of axial asymmetry) occur. The characteristic nuclear shapes occuring 
in the IBA are depicted at the vertices of the symmetry triangle [1] of the 
model (Fig. 1), labeled by their underlying dynamical symmetries, which are 
i) U(5), corresponding to near-spherical (vibrational) nuclei, ii) SU(3), repre- 
senting axially symmetric prolate deformed (rotational) nuclei, and iii) 0(6), 
describing nuclei soft with respect to axial asymmetry (7-unstable). 

Shape/phase transitions from one nuclear shape to another were first discussed 
in the context of the IBA in Ref. [6j, applying catastrophe theory to the energy 
functional [7] obtained in the classical limit of the IBA Hamiltonian through 
the use of the coherent state formalism [8f9] . A first order phase transition (in 
the Ehrenfest classification) has been found between the limiting symmetries 
U(5) and SU(3), while a second order phase transition has been located be- 
tween U(5) and 0(6). The spherical and deformed phases are separated by a 
narrow shape coexistence region [7] (also shown in Fig. 1), shrinking into the 
point of second order phase transition as the U(5)-0(6) line is approached. 

More recently, shape/phase transitions have been considered also in the frame- 
work of the geometrical collective model, resulting in the introduction of the 
critical point symmetries E(5) [10] and X(5) [11]. E(5), which corresponds to 
the second order transition between U(5) and 0(6), is a special solution of the 
Bohr Hamiltonian ^ using a potential u(/3) independent of 7 and having the 
shape of an infinite square well potential in 13. X(5), which corresponds to the 
first order transition between U(5) and SU(3), is a special solution of the Bohr 
Hamiltonian using a potential of the form u{l3) +f (7), where u{j3) is the same 
as before, while ^(7) is a steep harmonic oscillator centered around 7 = 0. 
E(5) and X(5) are also shown in Fig. 1, close to the points of the second and 
first order phase transitions of the IBA, respectively. 



Table 1 

Order v for states with any J, and for the special case of J'^ = 0^ states, in the 
geometrical models E(5), X(5), Z(5), Z(4), and X(3). J is the spin of the level, 
r = J/2, and n^ is the wobbling quantum number [3j which is zero for 0"*" states. 
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Table 2 

(Left) Energies of 0"^ states in the E(5), Z(5), and X(5) models. Energies on the left 
are in units of E{21) = 1.0, while in the column Norm, in units E{^2) = 1.0. The 
normalized results are identical for each of the models. The column IBA-Norm gives 
the normalized 0"*" energies for a large Nb IBA calculation near the critical point. 
(Middle) Same for the Z(4) model. (Right) Same for the X(3) model. Adopted from 



Ref. 


5]. 
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The competition between regular behavior dictated by underlying symmetries 
and chaotic behavior corresponding to lack of symmetries, has been studied 
throughout the symmetry triangle of Fig. 1. A highly regular region has been 
found along the U(5)-0(6) line, expected because of the underlying 0(5) sym- 
metry known to be preserved along this line [12]. Quite surprisingly, another 
regular region [13], connecting U(5) to SU(3) and called the Alhassid-Whelan 
arc of regularity, has been found inside the triangle, as shown in Fig. 1. The 
symmetry underlying the arc is yet unknown. 

0+ states are particularly appropriate for the detection of underlying symme- 
tries, because of the lack of centrifugal effects. As seen in Table 1, the critical 
point symmetries E(5) and X(5) mentioned above, as well as the Z(5) model 
|14] [a solution of the Bohr Hamiltonian similar to X(5), using an infinite 
square well in /?, but with t>(7) centered around 7r/6] possess as eigenf unctions 
the Bessel functions Jv While for J 7^ the order u of the Bessel functions 
is different in each solution, for J = the same order is obtained in all three 



Table 3 

0"*" bandheads in various (A, fi) irreps of the SU(3) limit of IBA. N stands for the 

boson number, Nb- Adopted from Ref. [5j. 
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(2N-12,6) (6N-15)/(2N-1) (2N-10,2) (6N-10)/(2N-1) 

(2N-16,8) (8N-28)/(2N-l) (2N-14,4) (8N-21)/(2N-1) 

Table 4 

0^ bandheads in various (o") irreps of the 0(6) limit of IBA. A^ stands for the boson 

number, Nb- Adopted from Ref. [5]. 



irrep 0"^ irrep 0^ irrep 0"*" irrep 0^ irrep 0^ 

(N) (N-2) 1 (N-4) 2 (N-6) 3-3/N (N-8) 4-8/N 

cases. As a result, 0^ states in these models look different if normalized to the 
energy of the 2f state, but they become exactly identical if normalized to 02", 
as seen in Table 2. Going further, one sees that in the latter normalization the 
energies of 0^ states, where n is their ordinal number, follow the simple rule 
n{n + 3). This is due to the fact that the spectrum of the roots of the Bessel 
functions J^, follows the n{n + z/ + 3/2) rule to a very good approximation for 
low u, being exact for z/ = 1/2 [5J. As a consequence, 0^ states in the Z(4) 
model [15] [similar to Z(5), but with 7 fixed to tt/Q] follow the rule ?t,(?t, + 2.5). 
Also, 0,^ states in the X(3) model [16] [similar to X(5), but with 7 fixed to 0] 
follow the rule n{n + 2), as seen from Tables 1 and 2. 

Taking into account the second order Casimir operator of the E(5) algebra 
[TSfTT] . the Euclidean algebra in 5 dimensions, one can see that the 0^ states 
in X(5) and Z(5) represent a case of a partial dynamical symmetry [IB] of Type 
I [in], a situation in which part of the states (the O"*" states in the present case) 



preserve the whole symmetry. 

What is a nontrivial result [S] , is that an IBA calculation near the point of the 
first order phase transition leads to a spectrum of 0^ states also following the 
n(n + 3) rule, dictated by infinite well potentials used in the Bohr Hamiltonian 
utilizing 5 degrees of freedom (the collective variables /3, 7, as well as the three 
Euler angles), as seen in Table 2. IBA calculations are performed using the 
usual IBA Hamiltonian [2U], involving two parameters {(, x)- Large boson 
numbers can be reached using the recently developed IBAR code [2II22]- 



It should be noticed that the n{n + 3) behavior of the 0^ states of the IBA 
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Fig. 2. Energies of low-lying states (normalized to E{2f)=l) of the usual IBA Hamil- 
tonian [20] with x=— \/7/2, C=0.4729, and Nb=250. C was chosen to reproduce the 
approximate degeneracy of E{02) and E{6^). Adopted from Ref. |23j . 

near the critical point of the first order phase transition is very different from 
the behavior of 0^ states obtained in the three dynamical symmetries of the 
model. Indeed, in the U(5) limit the energies of O"*" states increase linearly with 
the number of d bosons, i.e. with the phonon number. Analytical results for 
O"*" bandheads in the SU(3) and 0(6) limits of the IBA are shown in Tables 3 
and 4. In both cases in the limit of large boson numbers a linear increase is 
obtained. We conclude that in all three dynamical symmetries of the IBA, 0^ 
bandheads in the large boson number limit increase linearly, E = An. 

The regular behavior of 0^ states near the point of first order phase transition 
in the IBA invites a search for regularities of states with nonzero angular 
momenta. Indeed, as seen in Fig. 2 in an IBA calculation near the critical 
point, the O"*" bandheads are approximately degenerate with alternate levels 
of the ground state band with J/2 odd. 

Further investigation [23] of these degeneracies shows that the locus of the 
degeneracy -E(O^) = E{6f) [which is a hallmark of X(5)] in the IBA symmetry 
triangle is a straight line approaching the coexistence region in the limit of 
large boson numbers, as seen in Fig. 3. Similar results are obtained for the 
degeneracies -£^(03") = E{10f) and -£^(04") = ii^(14^). One concludes that these 
degeneracies characterize the coexistence region until the U(5)-0(6) line is 
approached. 

The ratio E{6f)/E{02) [23j, related to the first of the degeneracies mentioned 
above, turns out to be a simple empirical order parameter able to distinguish 
between a first order and a second order phase transition. Indeed, as seen in 
Fig. 4, this ratio exhibits the same behavior as the order parameter 1/3 used in 
Ref. [23]. A first order phase transition is seen for x = ~"\/7/2, while a second 
order phase transition is seen for x = 0- The ratio exhibits a sharp maximum 
just before the critical point in the first case (the effect becoming stronger at 
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Fig. 3. (Left) Line of degeneracy between the O2 and QJ levels for Nb = 10, 40, 
100, and 250 in the IBA triangle. (Right) Line of degeneracy between the Oj^ and 
10^ levels for Nb = 250 (top) and between the 0| and 14^^ levels for Nb = 250 
(bottom) in the IBA triangle. The dashed lines denote the critical region in the 
IBA obtained in the large Nb limit from the intrinsic state formalism. Adopted 
from Ref. 1231. 
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Fig. 4. The ratio E{6f)/E{02) ^s a function of ( for three values of x for (a) Nb 
15 and (b) Nb = 100. The inset to (a) shows the corresponding behavior for 1^2 
Adopted from Ref. [23]. 



larger boson numbers), while in the second case its behavior is smooth. 



Experimental data around the A^ = 90 isotones, the best empirical examples 
of X(5) [25|26ll27f28f29j . do exhibit in Fig. 5(a) a clear maximum just before 
A^ = 90, in agreement with the behavior expected for a first order phase 
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Fig. 5. (a) Experimental E{^{) j Ei^^^ ratio as a function of neutron number for 
the Nd, Sm, Gd, and Dy isotopes, (b) Same for the Xe and Ba isotopes. For smaller 
neutron numbers, the O3 state was taken in the ratio if its B{E2) decay was con- 
sistent with the a = N — 2 state. This corresponds to A^ = 74 in Xe and N = 76,78 
in Ba. Valence (hole) neutron number increases to the left. Adopted from Ref. 



transition. In contrast, experimental data around ^^^Ba, the best example of 
E(5) [30)129] . shows in Fig. 5(b) the smooth behavior expected for a second 
order transition. 

In conclusion, the large boson number limit of IBA reveals many regularities 
for O"*" states, the ones near the critical point of the first order phase transition 
being very close to the behavior obtained in critical point symmetries in the 
framework of the Bohr Hamiltonian. Degeneracies of O"*" states with states of 
nonzero angular momentum turn out to characterize the coexistence region 
separating the spherical phase from the deformed one, while ratios of energies 
of such pairs of states, like E{6i)/E{02), turn out to serve as order parameters 
able to distinguish between first order and second order phase transitions. 
These degeneracies call for further investigations into finding the symmetries 
underlying them. The recent conjecture [HI] of a partial SU(3) dynamical 
symmetry underlying the Alhassid-Whelan arc of regularity is also receiving 
attention. 

Work supported in part U.S. DOE Grant No. DE-FG02-91ER-40609 and under 
Contract DE-AC02-06CH11357. 
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